Spatial properties of entangled photon pairs generated in nonlinear layered structures 
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A spatial quantum model of spontaneous parametric down-conversion in nonlinear layered struc- 
tures is developed expanding the interacting vectorial fields into monochromatic plane waves. A 
two-photon spectral amplitude depending on the signal- and idler-field frequencies and propagation 
directions is used to derive transverse profiles of the emitted fields as well as their spatial correla- 
tions. Intensity spatial profiles and their spatial correlations are mainly determined by the positions 
of transmission peaks formed in these structures with photonic bands. A method for geometry op- 
timization of the structures with respect to efficiency of the nonlinear process is suggested. Several 
structures composed of GaN/AIN layers are analyzed as typical examples. They allow the generation 
of photon pairs correlated in several emission directions. Photon-pair generation rates increasing 
better than the second power of the number of layers can be reached. Also structures efficiently 
generated photon pairs showing anti-bunching and anti-coalescence can be obtained. Three reasons 
for splitting the correlated area in photonic-band-gap structures are revealed: zig-zag movement of 
photons inside the structure, spatial symmetry and polarization-dependent properties. Also spectral 
splitting can be observed in these structures. 

PACS numbers: 42.65.Lm,42.70.Qs,42.50.Dv 



-I— > 



> 



I. INTRODUCTION 

The first temporal correlations between the signal and 
idler photons generated in the process of spontaneous 
parametric down-conversion (SPDC) were observed al- 
ready more than thirty years ago They manifested 
quantum entanglement in the common state of both pho- 
tons generated in one quantum event of the spontaneous 
nonlinear process Since then the understanding of 
properties of such photon pairs has grown enormously. 
Polarization properties and namely polarization entan- 
glement between the signal and idler photons originat- 
ing in tensorial character of the nonlinear susceptibility 
have attracted the greatest attention from the very be- 
ginning. The reason lies in the simplicity of such states 
described in the Hilbert space with dimension 2x2. This 
made these states easily experimentally accessible. De- 
spite their simplicity, such states have allowed to reveal 
many substantial features of quantum physics related to 
correlations between subsystems [3| . Violation of the Bell 
inequalities that ruled out neoclassical theories [4| , the ef- 
fect of collapse of a wave-function 4j and the ability of 
teleportation of a quantum state [5J belong to the most 
important. 

In the time domain, both photons occur in a very nar- 
row temporal window [l[ because they are emitted at one 
instant after the annihilation of a pump photon. More- 
over the conservation law of energy dictates entangle- 
ment among monochromatic components of the signal 
and idler fields 0j. This entanglement leads to a typical 
finite entanglement time between the detection instants 
of the signal and idler photons observable in a Hong-Ou- 
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Mandel interferometer MJor sum- frequency generation 
of the paired photons [g, 0] . These properties have been 
extensively studied namely in connection with the gen- 
eration by pulsed pump fields [H, |i| that allows precise 
synchronization of photons from different photon pairs. 

Spatial properties of photon pairs have attracted at- 
tention last. Correlations in the transverse planes of 
the signal and idler photons occur here due to geometric 
properties of a photon-pair source and pump-beam spa- 
tial profile. They originate in the necessity of reaching 
good spatial phase matching of the interacting fields [2| 
that results in an efficient nonlinear process. For exam- 
ple, the sum of wave vectors of the signal and idler fields 
has to approximately give the wave vector of the pump 
field to observe an efficient photon-pair generation in an 
extended bulk crystal using collimated pumping. This re- 
sults in strong correlations in emission directions of the 
signal and idler photons [Iol - [l2l ] . These correlations may 
even be exploited to 'transfer' spatial properties of the 
pump beam into spatial correlations of the signal and 
idler beams [H, [Tj]- We may observe analogy between 
spatial and spectral correlations (entanglement) of two 
photons in a pair. Similarly as spectral correlations may 
be tailored by a pump-field spectral profile, spatial cor- 
relations may be controlled by a pump-beam spatial pro- 
file. Also signal and idler electric-field phase variations 
in the transverse plane belong to important character- 
istics of photon pairs. They can be quantified in terms 
of eigenstates of the angular orbital momentum operator 
[15| . Under certain conditions, entanglement between 
such states in the signal and idler fields has been ob- 
served [l6], [I?} ■ Experimentally, spatial correlations can 
be conveniently measured by moving a fiber tip in the 
transverse plane, as it has been done, e.g., in [18] . The 
use of an intensified CCD camera has provided a more 
elegant way for the experimental investigations of spatial 
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correlations [13, l20|. Spatial correlations between the 
down-converted beams have also been extensively used 
for quantum ghost imaging 21]. We note that spatial 
correlations occur not only in spontaneous regime, they 
have been found also in the regime with prevailing stim- 
ulated emission [H, [HI • 

All these forms of entanglement are potentially inter- 
esting both for fundamental physical experiments and 
practical applications including metrology 24] , quantum 
cryptography [25[ and quantum-information processing 
(3|. In principle, all these forms of entanglement can 
occur simultaneously depending on the source of pho- 
ton pairs. However, usual sources of photon pairs are 
constructed such that only one or two forms of entangle- 
ment (e.g., polarization and spectral) are found and may 
efficiently be experimentally exploited. 

The effort to generate entangled photon pairs extended 
over as many basis states as possible (and defined in sev- 
eral degrees of freedom) and to modify 'the structure of 
entanglement' belong together with the effort to enhance 
quantum efficiency of the generation process to the leit- 
motifs of the development in this field. Bulk nonlinear 
crystals that were nearly exclusively used at the begin- 
ning have been gradually replaced by more complex and 
efficient nonlinear structures including poled nonlinear 
materials [f| [26l - [28| . nonlinear waveguides (29j - [3l| and 
nonlinear photonic structures [32| . Nonlinear photonic 
structures are extraordinarily interesting because they 
allow efficient photon-pair generation owing to enhanced 
electric-field amplitudes observed inside structures with 
photonic band-gaps [H, H3] on one side, they also allow 
relatively wide tailoring of prop erties of the emitted pho- 
ton pairs on the other side [3 5m 71 ] . Structured nonlinear 
fibers that rely on four-wave mixing [38-40] represent 
a ty pical example. Or waveguides with Bragg mirrors 
[4l[ |42j can be mentioned as perspective sources. 

In the article, we consider systems composed of paral- 
lel nonlinear layers. Back-scattering of the fields on the 
boundaries between layers provides suitable conditions 
for the enhancement of electric-field amplitudes under 
certain conditions [43| . Moreover, spatial properties of 
photon pairs can be efficiently tailored using parameters 
of these systems (e.g., the number of layers). Nonlinear 
layered structures have already been studied in the sim- 
plified geometry both in the framework of classical (36| 
and quantum [3 71 ] descriptions. Here, we generalize the 
quantum model present in [37j to the real spatial geome- 
try including vectorial character of the interacting fields. 
This allows us to study transverse intensity profiles of the 
down-converted beams as well as correlated areas of the 
signal and idler photons considering typical layered struc- 
tures made of GaN/AlN. We note that GaN/AIN struc- 
tures with random layers' lengths have already been stud- 
ied as sources of photon pairs with ultra-narrow spectral 
widths that originated in an optical analog of Anderson 
localization of the down-converted fields |4|| . It has 
also been shown that GaN/AIN nonlinear layers allow 
the generation of photon-pair states antisymmetric with 



respect to the exchange of the signal- and idler-field fre- 
quencies that exhibit anti-bunching and anti-coalescence 

The investigated layered structures as sources of pho- 
ton pairs can be compared with other photon-pair sources 
with respect to intensity transverse profile, correlated 
area and efficiency of photon-pair generation as follows. 
The vast majority of photon-pair sources including wave- 
guiding structures and nonlinear crystals with/withour 
poling are designed such that the signal (and similarly 
the idler) photon is generated into one compact and small 
emission area. In some sources, bulk crystals are cut in 
such a way that the spatial symmetry allows to generate 
photons around the whole cone surface [Hj]. In this case 
the generated state is entangled also in the wave vectors 
of the signal and idler photons. As the Hilbert space cor- 
responding to this degree of freedom has a greater num- 
ber of independent states we obtain a higher-dimensional 
entangled state. In case of layered structures, the signal 
(and idler) photon can even be generated around sev- 
eral concentric cone surfaces depending on complexity of 
the structure. This even enlarges the number of inde- 
pendent states in the transverse area making the gener- 
ated state suitable for 'parallel processing' of quantum 
information in the transverse plane. As for correlated 
areas, they are described by compact profiles (usually of 
an elliptical shape) in all common photon-pair sources 
with cw pumping. Splitting of correlated areas into sev- 
eral well-separated parts is a distinguished property of 
nonlinear layered structures. Even three different mech- 
anisms leading to this splitting exist in layered struc- 
tures: zig-zag movement of photons inside the structure, 
spatial symmetry and polarization-dependent properties. 
The last two mechanisms may in principle occur also in 
other photon-pair sources. However, enhancement of the 
fields' amplitudes inside a structure is needed to make 
them significant. Finally, periodically-poled nonlinear 
crystals provide the greatest photon-pair fluxes. On the 
other hand, wave-guiding structures including nonlinear 
planar wave-guides and nonlinear structured fibers have 
the greatest nonlinear conversion efficiencies due to the 
transverse confinement of the interacting fields. How- 
ever, they allow for only moderate pumping intensities 
because of a possible material damage and competing 
nonlinear processes. Nonlinear layered structures lie in 
the middle in this comparison. They provide greater non- 
linear conversion efficiencies than bulk crystals owing to 
enhanced electric-field amplitudes along the propagation 
direction. Compared to wave-guiding structures, their 
nonlinear conversion efficiencies are lower because the 
'fields' confinement' is only in one dimension. On the 
other hand, similar pumping intensities as those used for 
bulk crystals can be applied. 

The paper is organized as follows. A spatial vectorial 
quantum model of SPDC in layered media is present in 
Sec. II and provides quantities characterizing the emitted 
photon pairs. A systematic method for designing efficient 
nonlinear structures is described in Sec. III. Transverse 
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intensity profiles and correlated areas are discussed in 
Sec. IV using structures with different numbers of layers. 
Conclusions are drawn in Sec. V. 



II. SPATIAL QUANTUM MODEL OF 
SPONTANEOUS PARAMETRIC 
DOWN-CONVERSION 









-11/ / 

1 


k 




1. 




e 


n (l) 








d (2) 




Zo 


Zl 





Nonlinear Hamiltonian Hj nt characterizing SPDC in a 
nonlinear medium of volume V at time t can be written 
in the form |2|: 



e / dr 
IV 



H int (t) 

JV 

d(r): fEW(r ) t)E(-)(r,t)E,H(r,t) + h.. 



(1) 



In Eq. ([T]) , d stands for a third-order tensor of nonlinear 
coefficients and the symbol : means shorthand of the ten- 
sor d with respect to its three indices. A strong pump 
field is described by the positive- frequency part Ep (r, t) 
of its electric-field amplitude vector. A signal [idler] field 
at single-photon level is characterized by the negative- 
frequency part Eg '(r, t) [Eg \r,t)} of its electric-field 
operator amplitude. Symbol eo denotes permittivity of 
vacuum whereas h.c. stands for the Hermitian conjugate 
term. 

The positive-frequency amplitudes E^ (r, t) of the in- 
teracting fields (m = p, s, i) can be in general decomposed 
into plane waves with wave vectors k m and amplitudes 
F^+Vk V 



EL +) (r,0 



1 d 3 k m E^ h) (k m )exp(ik„ i r-ia; m t); 

(2) 

uj m is the frequency of field m determined in accordance 
with dispersion relations. 

Considering the incident and un-depleted pump field, 
we assume that its temporal spectrum £ p (uj p ) as well 
as spatial spectrum £ p T (k px , k py ) in the transverse plane 
are given. In this case, the decomposition of amplitude 
E p (r, t) in Eq. @ can be rewritten as: 



E(+)(r,i) = 



1 



tt/2 pTr/2 

sin{'d p )d'& p i dtjj p 

-tt/2 



-tt/2 



dujp £ p (ujp)£ p [k PtX (r2p), k p y (J~2p) 



x exp [ikp^ x (fl p )x + ik Pi y(il p )y} 
x exp(— iuipt); 



q-TE,TM 



(3) 



using the vector O p = (uj p , ■d p , ip p ) of 'spherical coordi- 
nates' ujp, dp and tp p . Speed of light in vacuum is denoted 
as c. Assuming for simplicity air around the structure 
the x and y components of wave vector k p in front of the 



FIG. 1: Scheme of the structure and used coordinate system. 
A plane wave with wave vector k propagates along the radial 
emission angle ■& and azimuthal emission angle ip. The radial 
emission angle ■d is measured in the plane of incidence from 
the +z axis. The azimuthal emission angle ip gives rotation in 
the xy plane beginning from the +y axis and rotating towards 
the —x axis. Vectors exE and eiM give polarization directions 
of TE and TM waves determined with respect to the plane of 
incidence. Symbols Zi for i — 0, . . . , N identify z positions of 
the boundaries being perpendicular to the z axis; n"' [d"'] 
means index of refraction [tensor of nonlinear coefficients] in 
an Ith layer. 



structure are given as: 



kp.x (yip) 



ojp sin(^p) sin(# p ) 
c 

uip cos('0p) sin(?9p) 



(4) 



The decomposition of pump-field amplitude E p + ' into TE 
and TM waves as given in Eq. ([3]) is done with respect 
to the plane of incidence of the wave with wave vector 
kp propagating through the layered structure (see Fig. Q] 
for a scheme of the structure). We note that projections 
of k vectors into the planes of boundaries are conserved 
through the structure. The pump-field amplitudes ~E P ^2 
introduced in Eq. ^ describe field evolution along the 
z axis where the field undergoes back-scattering at the 
boundaries. 

Considering a structure with N layers and boundaries 
perpendicular to the z axis localized at positions z n , n — 
0, . . . , N, the pump-field amplitude E p ^ can be written 
in the form (for more details, see [37|): 



E p +)(z,Op) = rect 



E 

a=F,B 



x e 



(0) 



s (l) 



a (n p )exp [tKl°2 z (n p )(z - z ) 

(z) ]T < } , Q (Op) 

a=F.B 



N 

y ] rect Zi 



1=1 

-pJ,a( n p) ex P 

+ rect 



a=F,B 



4^ exp [iK^H^iz - z N ) 
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a = TE, TM. 



(5) 



The function rect ZajZ6 (z) equals one for z a < z < zi, and 

is zero otherwise. Symbols e P p ta and Gp\, a denote po- 
larization vectors of a waves for forward- and backward- 
propagating fields with respect to the +z axis, respec- 
tively. 

The z component Kp}_ z (Q p ) of wave vector in an ith 
layer belonging to wave a with frequency lu p and propa- 
gating in direction ($ p ,ip p ) in front of the structure can 
be expressed as: 



,(0 



■ cos 



(6) 



where the sign + (— ) is appropriate for a forward- 
backward-) propagating wave. Index of refraction of the 

pump field in an ith layer is denoted as n p . As we as- 
sume that the layered structure is surrounded by air, we 
have n p ' = n p N+1 ^ = 1. The propagation angle in 
the ith layer is then derived from Snell's law: 



rip ** sm 



( 1 ?(0) )=n (0 sill(l? W )) l = l,..., N +l; (7) 



i) 



(0) 



p — u p . 

The coefficients A P l^ a (ft p ) and A P g ta (fl p ) introduced 
in Eq. ([5]) determine amplitudes of a waves with fre- 
quency uip propagating forward and backward, respec- 
tively, in the direction. These coefficients are 
derived from Fresnel's relations at the boundaries using, 
e.g., the transfer-matrix formalism [43[. Values of the 

coefficients A$>, a (fl p ) and A p ^\fl p ) for a = TE, TM 
characterize the pump fields incident on the structure 
from both sides and represent the boundary conditions. 
The transfer-matrix formalism has been elaborated in de- 
tail in [37j for layered structures. The coefficients A® of 
the fields with wave vectors lying in the plane of incidence 
given by an angle ip p are given by Eq. (23) of Ref. [37| . 

The signal and idler fields are at single-photon level 
in the spontaneous process and so they have to be 
treated using quantum theory ||. Their electric-field 
operator amplitudes in the layered structure can be de- 
composed into plane waves using Eq. ([I]) and subse- 
quently described in analogy to the pump field. Using 
the 'spherical coordinates' u> m , i? m and ip m and defin- 
ing O m = (u> m , $ m , ipm), m — s,i, we can express the 

positive-frequency part Em ^ of the electric-field opera- 
tor amplitude of mode m as: 

, 1 r l% r /2 

E^M) = ^3 J / sin(i? ro )dtf ro / d^ m 



-7T/2 



ui m du} m exp [ik miX (n m )x + ik m ^(Vl m )y] 

E 



a=TE,TM 



167r 3 eo 



(z,n m )exp(-iLJ m t); (8) 



h stands for the reduced Planck constant and E 
ir(+)t 



(-) 



EJb . The expression ^hoj m / (167r 3 e ) gives an electric- 
field amplitude per one photon with energy huj m prop- 
agating at speed c. The x and y components of wave 
vectors k m of the signal and idler fields (m = s, i) are 
defined outside the structure: 



^m,rc(^m) 



u> m sin(^ 

III 

) sin(i? 

in 

c 

uj m cos(*0 m ) sin(i? m ) 



(9) 



The operator amplitudes a m!Q ,(z, fl m ) introduced in 
Eq. @ can be derived in the considered layered structure 
in the form: 



^m,a(z,^m) — rCCt_ oo _ 2:0 (z) ^ ^ Qj m ^ a (£l m *) 

a=F,B 

x e£) !a (n m )exp [iK£lja m )(z - z Q ) 

N 

+ ^2rect Zl _ 1:Zl {z) a m a ,a( n m) 



/=1 



-F,B 



x e®, a exp \iK m \ z {Vt m ){z - zi-! 
+ rect ZjVi00 (z) Ca'f^) 



a=F,B 



e^ a + 1 )(n m )exp [i^+ 1 )(n ro )(« - 2JV )' 



a = TE, TM. 



(10) 



Here, the symbols eS Fj(1 and eS s , stand for polariza- 
tion vectors of a wave of field m propagating forward 
and backward, respectively. The annihilation operators 
(fi m ) are defined at the end of an ith layer for the a 
wave with frequency u> m of field m propagating along the 
direction (i? OT iV'm) either forward (a = F) or backward 
(o = B). 

The z component K m \^ z (£l m ) of wave vector in the ith 
layer is determined as: 

= ± m( ro) ro c 0s( ^)) 



m = s, 1; 



■cos^ m „ 
a = F,B. (11) 



The sign + (— ) stands for a forward- (backward-) prop- 
agating wave and n m gives index of refraction of field 
m in the ith layer. Snell's law provides the propagation 
angles -dm needed in the ith layer: 

nW( Wro )irin(0<fl) = 4f+ 1 )(c m )sin(^ +1) ), 

l = 0,...,N- (12) 



I'm — u r 



The operators aS t ,tt(^m) f° r i = 0, ...,AT+1 and 
fixed m, a, w m , i? m , and ip m are connected by unitary 
transformations at boundaries (Fresnel's relations) and 
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free-field propagation transformations. This means that 
the usual boson commutation relations obeyed by the in- 
cident fields are 'transferred' through the structure and 
the only nonzero commutation relations (47j are the fol- 
lowing ones: 

7 c 2 

[®m a ,a(^m)> ®m' a , ,a' i^m')] = | gm ^ )|o; 2 ^ m - m '^ a ' a ' 

S a , a 'Sij'S(uj m - Lo' m ,)5(d m - i)' m ,)S(ip m - tp' m ,). (13) 

The transfer-matrix formalism allows to express the op- 
erators am a ,a(£l m ) in terms of the operators ai^a (n m ) 
and ainj, j0 (n m ) describing the outgoing fields. The ap- 
propriate relations valid for fields with wave vectors lying 
in the plane of incidence given by an angle ip m are given 
in Eq. (21) of Ref. 

The generation of a photon pair in the process of SPDC 
governed by the Hamiltonian H mi in Eq. ([1]) is described 
by a perturbation solution of the Schrodinger equation 
to the first order at time t — > oo assuming an incident 
vacuum state |vac) in the signal and idler fields at time 
t — > — oo. The resulting state \ip)° n i can be derived in 
the form: 



\ out 



. N 

1=1 a,b,g=F,B q,^,7=TE,TM 



E 



n 



ir/2 



tt/2 



I: 



7T/2 



tt/2 



dljjr, 



iol-,dbJ n 



y/W s Ui £p(Up)£p(kp,x, k P:V )S(uJ p - UJ S - CJj) 
X S \kp^ x k s x ^i,a;] & \J^p,y ^s,y 



X 



d« : e W, a (n ro )eW* (n.)eW;(n 



x exp 
x L/sinc 



2 PaSb^g.Z 



st,,/3\--»;~i g ,7 

(n p , n s , f2i)i; 



i 



AK 

2 PaS|,l g ,Z 



(Jl p , f2 s , Cti)Li 



x 4i) >a (n ro )a^(n s )a« t 7 (^)|vac) 



(14) 



The functions Ai^^^, n s , Q<) = i<XU n P ) ~ 

iQ2z(^s) — stand for phase mismatches in 

the Ith layer. Symbol means the length of Zth layer 
(Li = z\ — zi-i). The transverse wave vectors k m<x (tl m ) 
and k my (ft m ) for m = p,s,i are defined in Eqs. (U]) 
and (j9|) and characterize the fields outside the struc- 
ture. We note that the approach based on the solution of 
Schrodinger equation does not take into account surface 
SPDC that generates additional photon pairs [HI |49[ . 

The phase-matching conditions in the transverse plane 
xy are described by two S functions in Eq. (fl"4")l . that de- 
termine the emission direction {'di,ipi) of an idler pho- 
ton provided that the signal-photon emission direction 
(zJs, "0s) is given and the pump- field is in the form of a 
plane wave propagating in direction ($ p , ip p ). Simple ge- 



ometric considerations provide the following formulas: 



4>i 



ipp + arctan 



uj s sin(i? s )sin(V'p - ip s ) 



arcsm 



oj p sin(i? p ) — uj s sin(i? s ) cos('0p — ips) 
uj p sin(t9 p ) 



LOi cos(ip p - Ipi) 

LO s COs(lp p - 1p s ) 

Wi cos(ip p - Ipi) 



sin(# s ) 



(15) 



If the pump beam is focused, the strict phase-matching 
conditions in the transverse plane formulated in Eq. f| 15[) 
are relaxed and we arrive at correlation areas with finite 
spreads. 

The expression in Eq. (|14[) for the state IV')""* can be 
rewritten into the form containing only the outgoing cre- 
ation operators a^j^g (fi m ) and %iV^ m ) form = s,i 
and 13 = TE, TM using the formulas in Eqs. (21) and (23) 
of Ref. [37j • The outgoing operators can be finally trans- 
formed into the polarization basis of the detector using 
a suitable unitary transformation. We assume that the 
detection plane is perpendicular to the field propagation 
direction (i?, ip) and its s-polarization (denoted as _L) is 
parallel to the horizontal plane xz; p-polarization (de- 
noted as ||) is defined by the orthogonality conditions. 
Assuming field m (m = s,i) at frequency uj m and propa- 
gating along the angles $ m and ip m the needed i9 m - and 
i/Vn-dependent unitary transformation can be written in 
the form: 



a {N+1) rn ) 
a (N+1) rn ) 



a {a) rn ) 

a {0) (n ) 



COs(Cm) Sm(Cm) 
-Sm(C m ) COs(C m ) 



j_(n TO ) 

|| (Om) 



COs(Cm) Sm(Cm) 
-sin(Cm) COs(Cm)_ 



a 



,x(n m ) 
(n ro ) 



(m{flm,1prn) = arCCOS 



COS 



,|| 



x sign(^ m ), 



J I + sin 2 (tp m ) tan 2 (9 m 
/i m .), m = s,i. (16) 



The function sign gives the sign of its argument. 
The newly introduced annihilation operators a mb , a (n m ), 
m = s,i, b = F, B, a — \\ , _L, describe the signal and idler 
fields in the polarization bases connected with the detec- 
tors. 

The terms of state 1^)°"* in Eq. (pT4]) describing the 
created photon pair can be decomposed into four groups 
according to the propagation directions of the signal and 
idler photons with respect to the +z axis (FF, FB, BF, 
BB). Inside these groups there occur four contributions 
that differ in signal- and idler-photon polarization direc- 
tions (|| ||, || _L, _L ||, -L-L). Each contribution can be 
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written in the form: 



KIM) = n r' 2 ^ m )dd m r /2 



-tt/2 



/>oo 

Jo 

x exp[-i(k°f + k° ut )r] exp[i(a; s + Wi )t]; 

a,& = F,B; at,/9 = ||, JL - (17) 

The wave vectors k°"* and k° ut characterize free-field 
evolution of the emitted signal and idler fields, respec- 
tively, outside the structure. The functions (f>™f (f2 s , flj) 
introduced in Eq. (|17[) characterize completely properties 
of the generated photon pair. It gives the probability 
amplitude of having an a-polarized signal photon at fre- 
quency uj s propagating along direction (i^V's) together 
with its /3-polarized idler photon at frequency Ui prop- 
agating along direction (t?^ , -0^ ) at the output ab of the 
structure. 

Intensity spatial and spectral properties of a photon 
pair [37l ] can be conveniently derived from a density n°jf 
of the mean photon-pair numbers belonging to the state 
The density n"f is defined along the expression 



\n Sa , a mn ib A^iMtf,i b - (18) 



The photon-number density operator n rna ^ a {Tl m ) is ex- 
pressed as 



(l~2 m ) — a (S~2 m )tt maja (S7 r 



(19) 



Using Eq. (TTT)) the formula for density n"f of mean 
photon-pair numbers in Eq. (|18|) attains a simple form: 



n^(n.,n < ) = |^(n a ,n < )| 



(20) 



A density n"^ b (Q s ) of mean signal-photon numbers 
can easily be determined from the 



in the state 



density n" b defined in Eq. (|18p and using the relation in 
Eq. (TJOD: 



/tt/2 ptt/2 
-tt/2 



tt/2 Jo 



(21) 



If spectral resolution in detection of the signal-field 
transverse profile is not available, a spatial density 
n* r '^'^(i? s , V>s) of mean signal-photon numbers emitted 
in the direction ( , d s ,ijj s ) is a useful characteristic. It can 



be derived using the density n"'^ b of mean signal-photon 
numbers given in Eq. (|2 1 1) : 

POO 

v>.) - / ^<f b (a). (22) 

Jo 

The above densities can be analogously defined also for 
the idler field. 



Spatial correlations between the signal- and idler-field 
photon numbers in their transverse planes can be quan- 
tified in terms of the fourth-order correlation functions 
n c ° b ' a '^ ("dsjips, fli, V'i) that give joint densities of photon- 
pair numbers such that a signal photon propagates along 
the ($ s , tp s ) direction and the idler twin propagates along 
the {di^tpi) direction: 

poo pOO 

JO JQ 

(23) 

Provided that the signal-photon propagation direction 
(tfsiV's) * s §i ven > the joint density n c ° b ' a '^ (-d sl ip®, ipi) 
remains a function of the idler-field emission angles 
and ipi and its profile defines a correlated area. The cor- 
related area of an idler photon determines an area in the 
idler-field transverse plane where an idler photon can be 
expected provided that its signal twin has been detected 
in the ($°,?/>°) direction. 

Finally, an overall mean photon-pair number N" b ^ re- 
lated to the state l^)"'^. can t> e found using the following 
relation: 



N 



a,0 
ab 



/tt/2 /.717 
sm{$ m )d$ m \ 



/2 poo 

dtprn / diUr, 
/2 JO 



n°f(n s ,ni). 



(24) 



An important feature of layered structures is an in- 
crease of the efficiency of nonlinear process due to the 
enhanced electric-field amplitudes caused by interference 
along the z axis. This increase can be quantified with re- 
spect to a certain reference structure which fully exploits 
the nonlinearity, but does not rely on interference. In 
this reference structure, there occurs no back-scattering 
of the propagating fields and also the nonlinear process 
is assumed to be fully phase matched. The orientations 
of nonlinear layers and polarizations of the interacting 
fields are such that the most intense nonlinear effect oc- 
curs. This reference structure provides an ideal photon 
pair with a signal photon emitted in an arbitrary direc- 
tion ( r & s ,ip s ) and an idler photon in the corresponding 
direction. Fixing -d s and ?/> s , the emitted pair can 
be described by the following output state \ipY s e l [com- 
pare Eq. (fli])]: 



= / u s duj s j ^duji^/uJsUJ 



N 



£ p (uj s + u)j) 2J max(d ( ' ) )Lid\ {u s )a\ (uj) | vac) , (25) 



1=1 



where dl(u> s ) [d\(uji)} stands for a signal- [idler-] field 
creation operator outside the reference structure. The 
function max used in Eq. (|25|) gives the maximum value 
among the elements of tensor d^. 

Using the reference structure a relative density 
Vs'ab(^ s ) °f mean signal-photon numbers belonging to 
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the state 



can naturally be denned as 



(26) 



where the signal-field photon-number density w"„ 6 is 
written in Eq. (|21l) . The reference signal- field photon- 
number density n r s cl characterizes the state iV^g 6 ! m 
Eq. (f2"S")) and does not depend on the propagation angles 
i? s and ip a . 

Properties of photon pairs in the time domain are 
complementary to those found in the spectral domain 
and belong to important characteristics of photon pairs. 
We can mention signal- and idler-field photon fluxes or 
coincidence-count interference patterns in different kinds 
of interferometers as examples. These photon-pair prop- 
erties can be investigated, e.g., using the formulas con- 
tained in Sec. IIC of Ref. [37| even in this spatial vectorial 
model. 

Numerical calculations that follow are performed for 
a cw pump field with a Gaussian transverse profile de- 
scribed as follows: 



f tr (k h 



/2vr 



exp 



kl) 



(27) 



In Eq. (|27|l , £ p gives an amplitude of the pump field with 
the carrying frequency lu®] r p denotes the width of am- 
plitude transverse profile. Normalization of the function 
£p is such that f dk x J dk y \£^{k x , k y )\ 2 = 1. In case of 
cw pumping, there occur formulas in the above equations 
that contain a formal expression 6 2 (lu). This expression 
has to be replaced by the expression 2T/(2n)5(u), where 
the detection interval extends over (— T, T) and the limit 
T — > oo may be considered. 



III. DESIGN OF AN EFFICIENT LAYERED 
STRUCTURE 

Here we consider structures with odd numbers N of 
layers made of two kinds of materials. Layers of mate- 
rial b of length If, are sandwiched by layers of material 
a having lengths l a . Also pumping at a defined carrying 
frequency uj® and impinging on the structure at normal 
incidence is assumed. The down-converted signal and 
idler fields are assumed to have nearly degenerate fre- 
quencies. 

We suggest a method for designing an efficient layered 
structure from the point of view of, in general, three- 
mode nonlinear interaction. It is based on two observa- 
tions: 

• An efficient nonlinear process occurs provided that 
all three nonlinearly interacting fields lie inside 
their transmission peaks. This behavior originates 
in the fact that the electric-field amplitudes of 



monochromatic fields with frequencies in transmis- 
sion peaks are enhanced inside the structure owing 
to constructive interference of back-scattered light. 
It follows from the band-gap theory that the closer 
the transmission peak to a band gap, the greater 
the electric-field amplitudes. 

• As numerical calculations have revealed, the over- 
lap integral over the amplitudes of three interacting 
fields giving the strength of the effective nonlinear 
interaction [see Eq. (JTJ] vanishes if the signal- and 
idler-field amplitudes along the structure are the 
same. This means that photon-pair states degener- 
ate in frequencies, emitted in symmetric directions 
and having the same polarizations cannot be gen- 
erated. 

We note that these facts have been found crucial in de- 
signing layered structures efficient for second-harmonic 
generation 50]. Considering a collinear interaction, the 
requirements are even more strict because of only one 
propagation direction. This requires specific approaches 
that rely, e.g., on tuning the frequencies of transmission 
peaks by changing the index of refraction of one type 
of the layers [50]. Alternatively, non-collinear second- 
harmonic generation has been considered. In this case, 
there exists one free parameter (radial emission angle) 
that can be varied in order to find an efficient structure. 
This represents an equivalent problem to that considered 
here and can be treated by the developed systematic ap- 
proach. 

Returning back to the considered layered structures 
they are characterized by three parameters: number TV 
of layers and lengths l a and lb of these layers. A detailed 
inspection has revealed that the number N of layers sig- 
nificantly determines the number of generated photon 
pairs as well as angular extensions of the emission ar- 
eas of photons in a pair (see Fig. [3] bellow). The greater 
the number N of layers the greater the number of gen- 
erated photon pairs and also the smaller the angular ex- 
tensions of emission areas. From practical point of view, 
the number N of layers is approximately fixed consider- 
ing these dependencies. Having the number N of layers 
fixed, there remain two adjustable parameters - layers' 
lengths l a and If,. However, these two lengths cannot be 
chosen arbitrarily because the pump field at the carrying 
frequency to® has to be in a transmission peak. 

Let us fix the lengths l a and h (together with the num- 
ber iV of layers) for a moment and determine the spec- 
tral intensity transmission T(lu) along the +z direction 
(of the pump- field propagation) using, e.g., the transfer- 
matrix formalism [431 ] . Increasing the frequency u> there 
occur forbidden bands one following the other. The dif- 
ference in central frequencies of the neighbor forbidden 
bands is roughly the same in accord with the band-gap 
theory. This may be convenient for spectrally nearly de- 
generate SPDC provided that the process of SPDC can 
be tuned such that the pump-field frequency lies inside 
a transmission peak close to the second forbidden band 
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whereas the signal- and idler-field frequencies are inside 
the transmission peaks near to the first forbidden band. 
We note that all three interacting fields can accommo- 
date themselves either to the transmission peaks above 
or bellow the forbidden bands. This helps to obey the 
quite strong requirements of the three-field interaction 
by varying the signal-field emission angle However, 
suitable conditions can only be revealed numerically. 

Now we return back to the problem in which the pump- 
field carrying frequency cu® is fixed. As follows from the 
considerations of the previous paragraph, useful struc- 
tures are those that have the first upper or the first lower 
transmission peak near to the second forbidden band at 
the pump- field carrying frequency w". Inspection of the 
behavior of structures with different layers' lengths l a and 
lb has shown that there exists a system of curves in the 
plane spanned by layers' lengths l a and lb that provides 
the required transmission peaks at the frequency D®. One 
curve corresponds to the first lower transmission peak of 
the second forbidden band. Similarly, another curve is 
associated with the first upper transmission peak of the 
second forbidden band. They can be revealed as follows 
using the scaling property of diffraction phenomena in 
optics. 

We consider a fictitious dispersion-free structure with 
the indexes of refraction appropriate for the pump-field 
frequency id® and define the corresponding optical lengths 
l° pt and ^ pt . It can be shown that the ratio L = 
^ pt /?° pt of optical lengths represents a suitable vari- 
able for parametrization of these curves. Suitable op- 
tical lengths ^° pt and l^ pt for a given value of the ratio 
L can be revealed easily using the scaling property of 
diffraction phenomena. We fix the value of length l° pt 
to, e.g., l° pt '° = A°/2 = ttc/w°. This gives the unit 
length of diffraction phenomena that can be, in princi- 
ple, chosen arbitrarily. As the ratio L is given, the opti- 
cal length /° pt '° is derived as ;° pt, ° = Ll° pt >° and we can 
calculate the spectral intensity transmission T p (oj p ) for 
this structure. We further identify the frequency w™ ax 
of the first lower (or upper) transmission peak of the 
second forbidden band. Then we have to 'transfer' the 
actual frequency w™ ax of the transmission peak to the 
required frequency to® using the scaling property. The 
scaling property provides the layers' optical lengths in 
the form: ^° pt = Z° pt ^ w°/^ nax and Z° pt = Ll opt . We 
note that the obtained lengths differ for the lower and 
the upper transmission peaks. 

In the next step we move along the obtained two 
curves (for the upper and the lower transmission peaks) 
parameterized by the ratio L and numerically analyze 
the structures. The maximum ry™ ax of relative density 
rj s (tu s , i9 s ,ip®) of mean signal-photon numbers taken over 
the signal-field frequency lo s and radial emission angle 
•d s assuming the fixed signal-field azimuthal emission an- 
gle 4>° s has been found a suitable quantity for monitoring 
efficiency of the nonlinear process. The signal- and idler- 
field polarizations are assumed to be fixed. The greater 
the value of maximum ry™ ax of relative density the closer 



the signal- and idler-field transmission peaks to the first 
forbidden band. A curve giving the dependence of max- 
imum ?7™ ax of relative density on the ratio L is thus a 
good indicator for choosing suitable layers' lengths. We 
note that the curve depends on polarization properties 
of the signal and idler fields as well as the azimuthal 
signal-field emission angle tp®. It is also possible to mon- 
itor another quantity in this procedure, e.g., the over- 
all mean photon-pair number N given in Eq. (1241) . Or 
the density of mean signal-photon numbers emitted for 
a fixed azimuthal emission angle "0° and determined as 

sin(z9 s )dz9 ;s n* r (i? s , -0s) using the formula in Eq. (12"2")) 
may be considered. 

As an example, the dependence of maximum ■qj-^ ,ma ' x 
of relative density of two photons propagating along the 
+z axis on the ratio L for the first lower and the first up- 
per pump-field transmission peak of the second forbidden 
band of the structure composed of 11 and 101 layers, re- 
spectively, is plotted in Fig. [5] (for details, see below). 
The curve in Fig. |S{a) appropriate for the structure with 
11 layers is continuous and rather flat. This means that 
there exists a whole continuous set of structures giving an 
efficient nonlinear process. On the other hand, efficient 
nonlinear structures with N = 101 layers are found only 
in peaks of the curve in Fig. [U(b) . This is a consequence of 
complex interference of back-scattered light inside these 
structures. We note that the appropriate transmission 
peaks in all three interacting fields do not necessarily ex- 
ist for all values of the ratio L. Also the greatest values 
of maximum r/'t'^' max of the relative density are found 
for the values of ratio L in a certain restricted region. 
As the graph in Fig. [2jb) shows, the greatest values of 
maximum ^'" ,max occur around L = 0.5. 

The curves in Fig. [2] demonstrate a 'weak' increase of 
values of the maximum 7^ L, H ,max f relative signal-field 
density with the increasing values of the number N of lay- 
ers. As the densities n T s of the mean signal-photon num- 
bers of the reference structure are linearly proportional 
to the second power of the length of nonlinear material 
inside the structure, the densities n s of the mean signal- 
photon numbers increase better than N 2 . This is very 
important, because the mean number of emitted photon 
pairs increases less than the second power of the length 
of nonlinear material in other structures producing pho- 
ton pairs like bulk crystals, periodically-poled crystals or 
wave-guiding structures because of 'longitudinal' phase- 
matching conditions. Nearly linear dependence of the 
mean photon-pair numbers on the length of a nonlinear 
structure is commonly observed. 



IV. INTENSITY TRANSVERSE PROFILES 
AND CORRELATED AREAS 

We consider three examples of layered structures made 
of GaN/AIN that demonstrate typical features of the 
emitted photon pairs. They generate photon pairs with 
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FIG. 2: Maximum rj s of the relative density of mean 

signal-photon numbers for states with photons propagating 
forward as a function of the ratio L of optical layers' lengths 
for the first lower [upper] pump-field transmission peak of 
the second forbidden band for the structure with N = 11 
[N = 101] layers; = deg. 



different polarizations of the signal and idler photons. 
They are designed in such a way that three different 
mechanisms of splitting the correlated area are observed. 
The pump field propagating at normal incidence at the 
carrying wavelength A° = 400 nm is assumed. The signal 
and idler fields are emitted around the degenerate wave- 
lengths = X° L = 800 nm in non-collinear geometry. 
The GaN and A1N layers are positioned such that their 
optical axes are perpendicular to the boundaries. The 
considered structures differ in the number N of layers 
(N = 11, 51 and 101). This results in different signal-field 
intensity transverse profiles plotted in Fig. [3] Whereas 
the shortest structure with TV = 11 layers has only one 
emission area, the longest one with TV = 101 layers al- 
ready forms five emission rings. We note that transverse 
profiles of the down-converted fields can also be efficiently 
tailored by modifying the pump-beam profile fl4j |. 

The analyzed structures have certain properties in 
common. Namely, they cannot generate proton pairs in 
nearly collinear geometries because of the symmetry that 
cancels the overlap integral. As for correlated areas, their 
azimuthal spreads Aijji depend on the width r p of trans- 
verse pump- field profile [H, [5l[ . The radial spreads 
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FIG. 3: Densities nf of mean signal-photon numbers depend- 
ing on signal- field radial (i? s ) and azimuthal (tp s ) emission 
angles for (a) TV = 11, (b) TV = 51 and (c) TV = 101 layers 
without polarization resolution for both photons propagating 
forward. The emitted signal field is projected onto a hemi- 
sphere, one quadrant of which is plotted. In the graphs, the 
radial emission angle i? s determines the distance from the 
origin, whereas the azimuthal emission angle ip s gives the 
rotation measured from the vertical direction. Behavior of 
the density n* r in the remaining three quadrants can be de- 
rived from symmetry. The density n/ is normalized such that 
/; /2 aii(i)#,J° *n';(t)„i) = (tt/180) 2 /4; r p -> do. 



depend also on the geometry of the layered structure; the 
greater the number TV of layers, the smaller the correlated 
area. In more detail, the considered structures behave as 
follows. 
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FIG. 4: Relative density r] a of mean signal-photon numbers 
depending on normalized signal-field frequency 2uj s /uip and 
radial emission angle # s for the structure with N = 11 lay- 
ers for both photons with arbitrary polarizations propagating 
forward; ip s = deg, r p — » oo. 



A. Structure with 11 layers 

The structure is composed of 6 nonlinear GaN lay- 
ers 90.14 nm long and 5 linear A1N layers 74.92 nm 
long. Despite a relatively small number of layers, pho- 
tonic band gaps are already formed though the bot- 
toms of forbidden bands have intensity transmissions 
around 0.3. For ip B = deg, an efficient SPDC occurs 
for polarizations (TE,TM,TE) and (TE,TE,TM) for the 
(pump, signal, idler) fields due to the nonlinear coefficient 
d(l, 1,3). The fields are tuned to the first lower trans- 
mission peaks of the first and second forbidden bands, 
so the efficiency of SPDC is nearly optimal. The signal 
and idler transmission peaks are spectrally broad and, 
as a consequence, we observe one large emission area in 
the signal-field transverse plane that extends from cca 
■d s = 20 deg to 60 deg [see Fig. E^a)]. The relative den- 
sity rj s of mean signal-photon numbers as plotted in Fig.0] 
reveals that the emitted photons have nearly degenerate 
frequencies and relatively broad spectra. The shape of 
relative density rj s as shown in Fig. 2] is caused by a dif- 
ferent radial dependence of intensity transmissions Xte 
and Ttm- In the area around ip s = ±90 deg there occurs 
no photon-pair generation because of geometric reasons 
(the propagation of signal and idler fields as TM waves is 
not supported) . Compared to one GaN layer 6 x 90.14 nm 
long, the structure gives cca 2 times greater relative den- 
sities 7] s . 

The correlated area of an idler photon for a collimated 
pump beam is composed of several 'islands' localized 
around different idler-field radial emission angles $j [see 
Fig. 03a)]. This breaking is caused by the zig-zag move- 
ment of two photons inside the structure after being gen- 
erated in one layer. The comparison of correlated areas 
for collimated and focused pump beams using, e.g., the 
graphs in Fig. (5Ja) and (b) reveals that the spread of 
correlated area along the azimuthal emission angle Vi 
depends strongly on the amount of focusing the pump 
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FIG. 5: Correlated area n1 OT {-&i, ipi) of an idler photon for (a) 
r v — 1 mm and (b) r p = 30 nm belonging to a signal photon 
propagating along direction i)® = 38 deg and $J = deg. 
Both photons with arbitrary polarizations propagate forward; 
0i = tf? + 5&i, 0? = -0°, Vi = i>i + H, i>° = -V>°- The 
normalization f^^d&i f^^dipinf'^ijtpi) = (7r/180) 2 is 
used. 

beam. On the other hand, the shape of correlated area 
in the radial emission angle i?j is blurred by a focused 
pump beam keeping the overall spread of correlated area 
roughly the same. 

B. Structure with 51 layers 

The second structure contains 26 nonlinear GaN layers 
106.87 nm long and 25 linear A1N layers 65.99 nm long. 
This structure has already well-formed forbidden bands. 
As Fig. [UJb) shows there occurs an efficient SPDC in 
two concentric rings in the transverse plane. The struc- 
ture is designed such that an efficient nonlinear interac- 
tion for fields' polarizations (TM,TM,TM) occurs along 
the direction ip a = deg using the nonlinear coefficient 
d(2, 2, 3). The pump field lies in the first upper transmis- 
sion peak of the second forbidden band and the signal and 
idler fields are in the first and second lower transmission 
peaks of the first forbidden band. The relative density 
■q s of mean signal-photon numbers as plotted in Fig. [6] 
demonstrates that the signal-field (and also idler-field) 
intensity spectra are composed of two symmetric peaks. 
There occurs no photon-pair generation for degenerate 
signal- and idler-field frequencies because of symmetry 
This indicates the generation of photon pairs in the state 
that is anti-symmetric with respect to the exchange of the 
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FIG. 6: Relative density rj s of mean signal-photon numbers as 
it depends on normalized signal-field frequency 2u) s /u)p and 
radial emission angle # s analyzing the structure with N — 51 
layers and both forward-propagating photons with arbitrary 
polarizations; ip B = deg, r p oo. 



signal- and idler-field frequencies. Anti-bunching of pho- 
tons in a pair and anti-coalescence (observed in a Hong- 
Ou-Mandel interferometer) are the distinguished features 
of this state (46{. The graph in Fig. G^b) reveals that an 
efficient SPDC occurs also for azimuthal emission angles 
ips around ±90 deg. In this case, the signal photon is 
generated as TE wave and the idler photon as TM wave 
or vice versa. The structure gives cca 50 times greater 
relative densities -q s in comparison with one GaN layer 
26 x 106.87 nm long. 

As a consequence of two-peak structure of the signal- 
and idler-field spectra there occurs splitting of the corre- 
lated area of an idler photon, as documented in Fig.JTJa). 
The overall correlated area is composed of two symmet- 
ric parts. Because this splitting arises from symmetry, it 
survives even for focused pump beams, as demonstrated 
in Fig. [Tib). 



C. Structure with 101 layers 

The last structure is composed of 51 nonlinear GaN 
layers 106.42 nm long and 50 linear A1N layers 65.71 nm 
long. The pump field lies in the first upper transmission 
peak of the second forbidden band. A detailed compari- 
son of the relative density rj s of mean signal-photon num- 
bers (see Fig. [5]) with intensity transmission spectra Xte 
and Ttm valid for TE and TM waves, respectively, and 
obtained for i/j = deg (see Fig. [5]) reveals the follow- 
ing. An efficient SPDC occurs at intersections of jth and 
(j + l)th lower transmissions peaks of the first forbidden 
band for j = 2,3,4,5,6. This results in five concentric 
rings in the density n* r of mean signal-photon numbers 
clearly visible in Fig. (He). We note that the nonlinear 
coefficient (2(1,1,3) is exploited here and (TE,TM,TE) 
and (TE,TE,TM) [(TM,TM,TM)] polarizations are suit- 
able for the angles around tp s — deg [tjj s = ±90 deg]. 
The relative density r) s of mean signal-photon numbers 
indicates that the signal-field intensity spectra are com- 
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FIG. 7: Correlated area n™ r (i?i, ipi) of an idler photon for (a) 
r v = 1 mm and (b) r p — 30 nm corresponding to a signal pho- 
ton propagating along direction $® = 29 deg and ips=0 deg; 
N = 51. Both photons of arbitrary polarizations propagate 



forward; 0< = !?° + <5i?i, 




FIG. 8: Relative density r) B of mean signal-photon numbers 
depending on normalized signal-field frequency 2uj s /ujp and 
radial emission angle i? s . The structure with TV = 101 lay- 
ers and both photons of arbitrary polarizations propagating 
forward are assumed; ip s = deg, r p — s> oo. 



posed of two peaks of different weights for ijj s = deg. 
One peak is related to a TE signal-field wave, whereas 
the second one arises from a TM signal- field wave. Com- 
parison with one GaN layer 51 x 106.42 nm long reveals 
that the enhancement of optical fields inside the struc- 
ture results in an increase of the relative density rj s cca 
330 times. 

A two-peak spectral structure leads to splitting of the 
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FIG. 9: Intensity transmission coefficients Ttb for TE waves 
[(a)] and Ttm for TM waves [(b)] as they depend on radial 
emission angle i? assuming ip = deg and the structure with 
N = 101 layers. For ip = ±90 deg, the graphs for TE and 
TM polarizations are mutually exchanged. 
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FIG. 10: Correlated area n^^i, if)i) of an idler photon for 
(a) r p = 1 mm and (b) r v = 30 nm belonging to a signal pho- 
ton propagating along direction -d® = 23 deg and tp® = deg. 
Both photons have arbitrary polarizations and propagate for- 
ward; ■di = 0? + 5#i, = -0°, j>i = V? + Sipi, $ = -ip°. 



correlated area of the idler photon into two parts [52} , as 
shown in Fig. 110( a). However, two parts are not sym- 
metric in this case and they can also merge together 
provided that the pump beam is sufficiently focused [see 
Fig. fTCT b)]. The distance between two parts in the idler- 
photon radial emission angle i9j increases with the in- 
creasing signal-photon radial emission angle tf® . Whereas 
this distance equals cca 0.6 deg for fl® = 23 deg (the first 
ring), it equals already cca 4 deg for $° = 66 deg (the 
fifth ring) . Idler photons found in different parts of the 
correlated area differ in their polarizations. 

The above investigated structures demonstrate main 
features of photon pairs generated in nonlinear layered 
structures. These structures allow to generate photon 
pairs entangled in frequencies, polarizations as well as 
in emission directions. Moreover, they allow tailoring of 
properties of photon pairs varying namely the number of 
layers. A small amount of nonlinear material inside them 
is compensated by an increase of electric-field amplitudes 
originating in fields' back-scattering. 

V. CONCLUSIONS 

We have developed a spatial vectorial quantum model 
of spontaneous parametric down-conversion in nonlin- 



ear layered structures. Photon pairs generated in these 
structures can be entangled in frequencies, polarizations 
as well as emission directions. Namely entanglement in 
emission directions is important because it can be effi- 
ciently tailored varying the number of layers. A method 
for designing efficient layered structures has been sug- 
gested. Its efficiency has been demonstrated suggesting 
three typical structures. It has been shown that the num- 
ber of generated photon-pairs increases greater than the 
second power of the number of layers. Signal-field inten- 
sity profiles and correlated areas in the transverse plane 
have been analyzed in the designed structures. Inten- 
sity profiles are typically composed of concentric rings. 
The greater the number of layers, the greater the num- 
ber of rings. Correlated areas may be broken into sev- 
eral parts because of three possible reasons: i) Zig-zag 
movement of photons inside the structure, ii) Neces- 
sity to obey geometric symmetries, and iii) Polarization- 
dependent properties. Splitting of the correlated area 
arising from the geometric symmetry survives even for a 
focused pump beam. Also states of photon pairs showing 
anti-bunching and anti-coalescence can efficiently be gen- 
erated in layered structures. We believe that nonlinear 
layered structures are potentially interesting namely as 
efficient sources of photon-pair fields entangled in propa- 
gation directions which might be useful, e.g., in 'parallel 
processing' of quantum information. 
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